Collective Excitations of Trapped Imbalanced Fermion Gases 
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We present a theoretical study of the collective excitations of a trapped imbalanced fermion gas 
at unitarity, when the system consists of a superfluid core and a normal outer shell. We formulate 
the relevant boundary conditions and treat the normal shell both hydrodynamically and collision- 
lessly. For an isotropic trap, we calculate the mode frequencies as a function of trap polarization. 
Out-of-phase modes with frequencies below the trapping frequency are obtained for the case of a 
hydrodynamic normal shell. For the collisionless case, we calculate the monopole mode frequencies, 
and find that all but the lowest mode may be damped. 

PACS numbers: 



The recent pioneering experiments on polarized 
fermion systems in the BEC-BCS crossover make a broad 
range of exciting new phenomena accessible and offer the 
opportunity to study longstanding issues in the field of 
strongly interacting many-body systems [l], [1] ■ By load- 
ing a trap with unequal numbers of particles of two spin 
states, a phase-separated system consisting of an unpo- 
larized superfluid (SF) core and a surrounding polarized 
normal (N) phase is obtained. Experimentally probing 
the collective excitation frequencies in ultracold gases is, 
by now, a standard technique, and has proved indispens- 
able to understanding such systems. An experimental 
study of collective modes in imbalanced fermionic su- 
perfluid systems in the near future seems, therefore, a 
realistic prospect. 

The N-SF interface plays a major role in the static 
properties of the polarized system: taking into account 
the interface tension has proved essential to explain- 
ing the experiments performed in highly anisotropic 
traps [1, H]]. Moreover, its presence affects the ther- 
mal equilibration process at low temperatures Q and has 
been speculated to cause large fluctuations in the polar- 
ization at the interface [6[. 

In this Rapid Communication, we study the effects of 
polarization on the frequencies of the collective excita- 
tions. To do this, it is necessary first to establish the 
boundary conditions encoding the important physics at 
the interface. We then study two distinct cases: in one, 
the normal part is assumed to behave hydrodynamically, 
and in the other, collisionlessly, so that it is described 
by a Boltzmann-Vlasov equation. We present results for 
isotropic traps only; results for a highly elongated trap 
will be presented separately, along with more detailed 
calculations 

For a fully hydrodynamic trapped system, we find two 
general classes of excitations: A first class consists of 
modes which for zero polarization reduce to the well- 
known single-component-system excitations; for finite 
polarization, these evolve to modes with the same general 
character of motion but slightly shifted frequencies. In 
what follows, we refer to them as in-phase (IP) modes, for 
reasons that will become clear. Additionally, excitations 
unique to the two-component system constitute a second 



class. Particularly interesting is that, at low polariza- 
tion, these new modes have frequencies below the trap- 
ping frequency. Each of these low-frequency excitations 
is in one-to-one correspondence with the nodeless (n = 0, 
£) modes of the single-component system, to which their 
motion is related as follows: the outer boundary of the 
N and SF components oscillate in the same way as the 
corresponding one-component n = mode, except that 
they are out of phase with each other. This is illustrated 
in Fig. [1] for the quadrupole mode. For brevity, we re- 
fer to these new type of excitations as the out-of-phase 
(OOP) modes. 

In the collisionless case, we study the monopole (£ = 0) 
modes. Simultaneously solving the hydrodynamic and 
Boltzmann-Vlasov equations for the SF and N phase re- 
spectively, and imposing the appropriate boundary con- 
ditions, we find that the lowest, n = 1, monopole mode 
remains at 2u>q for all polarizations (here luq is the trap- 
ping frequency), in agreement with an exact result ob- 
tained by Castin [l8|. On varying the polarization, the 
n = 2 monopole mode frequency crosses over from the 
value appropriate for a fully-superfluid at zero polar- 
ization and a fully-collisionless system at complete po- 
larization. For this (and higher) monopole modes, and 
for some intermediate polarizations, damping occurs, the 
cause of which is purely geometric. 

We shall first describe the equilibrium state of the sys- 
tem, and then separately present our calculations and 
results for the hydrodynamic and collisionless cases. 

Equilibrium — Consider a spin mixture of polarization 
(JV T - N l )/(N^ + with N-\ and N t the spin up and 

FIG. 1: Time evolution of the quadrupole (I — 2) out-of-phase 
mode of an imbalanced fermion system in the hydrodynamic- 
normal-gas case". The outer boundaries of the normal (N) 
and the superfluid (SF) phase perform the same motion but 
7r out of phase. 
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down particle numbers, trapped by a harmonic poten- 
tial V(r). A nonzero imbalance generally results in the 
appearance of at least two phases: the unpolarized SF 
and the partially or completely polarized N phase which 
consists, in general, of a mixture of j and f particles. 

When the N and SF phases are separated by a first- 
order phase transition, an interface is formed between 
them with the SF in the trap center. Due to the uni- 
versal nature of the strongly interacting regime, the 
equilibrium pressure P and density p in the SF are 
known [9] to be related by P SF = fi 2 £(37r 2 ) 2 / 3 p^/5m, 
with the universal constant £ w 0.41 and m the parti- 
cle mass; in the N phase, on the other hand, one has 
Pi = fi 2 (67r 2 ) 2 / 3 p^ /3 /5rn for i =|, j (the overline de- 
notes equilibrium values). Moreover, within the local 
density approximation all three of f, I and SF may be 
described by an effective equilibrium chemical potential 
~pj(r) — (Xj — V(r), where is the chemical potential at 
the center of the trap and j =1,1 or SF (20j . 

Next we include the interface into our equilibrium 
framework. At zero temperature, the SF is theoreti- 
cally predicted to be unpolarized; experimentally, the 
trap core is observed to be unpolarized. This implies 
a negligible population of single-particle excitations in 
the trap center at experimental temperatures. Denoting 
the equilibrium position of the interface by £, chemical 
equilibrium is then ensured by imposing [la ] 

Furthermore, for mechanical stability, the Laplace con- 
dition must be satisfied: 




Here, i?i and i?2 are the radii of curvature of the inter- 
face, the position of which is denoted by £, and a are the 
interface tension 0, 0]; following Haque and Stoof [H, 
we take a = 0.6mfi'g F /H 2 [2l[. Note that psf is the to- 
tal chemical potential, including its fluctuation during an 
oscillation, which reduces to JigF at equilibrium. Equa- 
tions UJ and together with the expressions for the 
densities and pressures, allow us to fix the position of 
the interface for given particle number and polarization 
in the global equilibrium state. To do this, all that is 
necessary is to solve Eqs. (JlJ and @ for the radii while 
imposing the total particle number to be 7 x 10 6 . We 
then find that, for all polarizations, the superfluid core is 
surrounded by a fully polarized normal shell of the ma- 
jority species. This is similar to what is observed in the 
Rice experiments 0, H[ • 

Dynamical phenomena, such as the collective excita- 
tions studied here, involve situations in which the local 
velocities are nonzero. For such cases, involving depar- 
tures from the global equilibrium state, we must extend 
Eqs. fl]) and ^ appropriately as well as supplement them 
with additional boundary conditions. 
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FIG. 2: Collective mode frequencies against polarization for 
a trap with a hydrodynamic normal shell and with 7 x 10 6 
particles. The dashed and full lines denote the out-of-phase 
and in-phase collective modes respectively. The frequencies of 
the in-phase modes are rather insensitive to the polarization, 
except for polarizations after the avoided crossing. At zero 
polarization, the IP modes reduce to the excitations of the 
single-component systems, labelled by the numbers n and £ 

Dynamics: Hydrodynamic Normal Gas — Consider 
first the case in which the normal phase behaves hydro- 
dynamically, that is, when inter-particle collisions are fre- 
quent enough for local thermodynamic equilibrium to be 
ensured everywhere during the oscillation. The chemi- 
cal potential and pressure remain well-defined quantities 
and therefore Eqs. (T_TJ) and @, with 71 and P replaced by 
~p + 5/i and P + SP, respectively, remain applicable. 

In the system under study, particles may pass over 
from the N to the SF phase and vice versa Q . This inter- 
conversion separately conserves the mass of each species. 
Denoting the position of the interface by C = C+^Ci where 
5C is the departure of the interface from the equilibrium 
position, the appropriate boundary conditions (one for 
each of i = \ , j) are, at each point on the interface, 

e c • (2wipi - v S fPsf) c = (2pi - psF)(d t S(, (3) 

in which is the unit vector perpendicular to the in- 
terface (and directed towards the N side) and v are the 
velocities. These boundary conditions are, essentially, a 
restatement of the continuity equation for the case of two 
coexisting phases of different densities with a movable, 
permeable interface separating them. 

Finally, the bulk dynamics for both the SF and the N 
are governed by the Euler and continuity equations. Ap- 
plying standard techniques [ill ]. w ^ linearize these and 
obtain the equation for the deviation of the chemical po- 
tential from its equilibrium value, Sp. Taking the tem- 
poral and angular parts of SpsF, S^t, dpi, and <5C to 
be e luJt and F ; m respectively, with Y™ the spherical har- 
monics, we obtain 8/jsf <x r l V (a + , a - , a , (r/RsF) 2 ) 
and Sm oc r*F(a+, or, 3/2,1- {r/Ri) 2 ) with i =T,j 
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and F the hypergeometric function. Here a — I + 3/2, 
2a ± = i + 2 ± [P + I + 4 + 3w 2 /tj2]i/2 and £ a posi _ 

tive integer or zero. For all three phases, we have also 
defined R? = where j =1,1 or SF. The lin- 

earized versions of Eqns. {J), ([2]) and (J3]), together with 
the linearized bulk solutions, now suffice to describe the 
dynamical behavior of the system in the case where the 
normal side is hydrodynamic [Tj. 

Results: Hydrodynamic Normal Gas — In Fig. [2] we 
show the dependence of the mode frequencies on the po- 
larization. We distinguish between the two kinds of ex- 
citations already mentioned: the in-phase (IP) modes, 
analogous to the excitations present in single-component 
systems (full lines), and the out-of-phase (OOP) modes, 
unique to two-component systems (dashed lines). Note 
that similar excitations appear in trapped boson-boson 
systems 0. At zero polarization, the IP mode frequen- 
cies reduce to the single-component frequencies u> 2 /u)q = 
£ + 2n[(2/3)(ra +1 + 1/2) + 1] [3 which are in almost 
exact agreement with the experiments (l6| . For all po- 
larizations, the frequency of the sloshing (or Kohn) mode 
remains exactly u> = ujq; this excitation corresponds to 
a rigid cloud motion and is therefore not affected by the 
bulk equation of state. The frequencies of the IP modes 
with n = and I > 1, on the other hand, vary slightly 
with polarization; this is a consequence of the nonzero 
interface tension. We note here that recent experiments 
at MIT reveal the existence of a large partially polarized 
shell resulting from the interactions in the N phase. The 
incorporation of such interactions is beyond the scope of 
this work. Nevertheless, our main result, the existence 
of the OOP modes, remains valid there since such modes 
are generic to the two-component system (Toj . 

Dynamics: Collisionless Normal Gas — At tempera- 
tures well below the Fermi energy, a normal Fermi gas is 
expected to behave collisionlessly, especially in the case 
when this phase is fully spin polarized. We therefore 
turn our attention to a collisionless normal gas, that is, 
one that is described by a Boltzmann-Vlasov equation, 
as used, for example, to study a (spatially mixed) Bose- 
Fermi system in Ref. [l9|. The issue of collective exci- 
tations in normal fermion gases has been addressed us- 
ing several methods which include the hydrodynamic ap- 
proximation, the method of averaging, the scaling ansatz 
method, the sum-rule approach, and the random-phase 
approximation 

The fully polarized normal gas is now described by a 
distribution function /(r, v, i), evolving according to the 
Boltzmann-Vlasov equation 



dtf- 



V v f = 0, 



(4) 



rather than the hydrodynamic equations. We shall again 
consider small deviations from the equilibrium func- 
tion, taking / to deviate from the equilibrium valu e by 
a non-isotropic deformation of the Fermi surface [13j |: 
/(r,p,t) = /o(r,p,f) +S(\v\ -v F )v(r,p,t), with f the 
Fermi function and vf the (position-dependent) Fermi 
velocity. The interesting question then arises of which 
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FIG. 3: Monopole frequencies against polarization for a trap 
with 7x 10 6 particles. The lowest mode frequencies are exactly 
at 2wo, both for the cases of a collisionless (full line) and hy- 
drodynamic (dashed line) normal shell. The second monopole 
frequencies deviate substantially between these case. The 
grey full line indicates the region of damping for the colli- 
sionless case and the damping rates are shown in the inset. 



boundary conditions apply at an interface between a hy- 
drodynamic and collisionless gas |16j. First, notice that, 
since we no longer have local thermodynamic equilib- 
rium, there is no analogue of Eq. {T]) for the deviations 
from the equilibrium configuration. On the other hand, 
the Laplace condition, Eq. @ remains valid with the fol- 
lowing modifications: a) since the normal gas is fully po- 
larized, only quantities pertaining to the majority species 
appear in it; b) the role of the pressure fluctuation of 
the majority species, <5Pf in the previous section, is now 
played by the radial component of the momentum flux 
tensor 8Iil r = to 4 / d 3 «^(/ - /o)/(2tt?i) 3 . 

Next, from Eq. we see that the radial component of 
the SF velocity must be equal to the radial velocity of the 
interface, that is, • v$f = dtSC, so that there is no flux 
through the interface. There are therefore only two possi- 
bilities for a particle incoming on the interface from the N 
side: Andreev reflection or specular reflection Q. Here, 
Andreev reflection is suppressed by the full spin polari- 
sation in the N phase, leaving specular reflection as the 
sole reflection mechanism. According to Bekharevich and 
Khalatnikov, specular reflection off a moving interface re- 
sults in the boundary condition x) — v {x) ~ 2iLU\SC 
where \ = cos <f> and </> the angle between r and v. At 
r = C the rotationally symmetric (£ = 0) solution satis- 
fying Eqn. (U]) and this boundary equation is 



KC,X,*) = XW<5C[cot(wT/ 2 ) + »] 



(5) 



where r is the time for a classical particle of velocity v 
to travel from the interface and back f22[. 

Results: Collisionless Normal Gas — In Fig. [3J we 
present the monopole frequencies and compare the re- 
sults of the collisionless (full lines) with those of the hy- 
drodynamic approach for the normal shell. For both 
approaches, the lowest monopole frequency is 2ojq for 
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all polarizations; this is in agreement with an exact re- 
sult derived by Castin [lj| and differs from the results 
of Ref. [13] ■ As expected for the collisionless case, the 
second monopole frequency is •v/32/3o'o at zero polar- 
ization and 4u>o at full polarization. We find, however, 
that the collisionless mode is damped for polarizations 
between 0.09 and 0.28 (grey full line). The origin of this 
damping is a resonance which occurs when r = 2tt/lo; 
that is, when there exist particles for which the time to 
travel from the interface and back equals the period of 
the collective motion. Mathematically, it corresponds to 
the pole in v (see Eqn. ©). The damping frequency is 
obtained by writing lj — > u> — iT with to and T real and 
positive, and analytically continuing the pressure tensor 
IT rr from negative to positive values of T @, [H| • As is 
shown in the inset of Fig. [3J the maximal damping for 
this number of particles is T — 0.05wo, which is weak but 
experimentally detectable. We stress that the origin of 
this damping is purely geometrical, and its presence and 
strength will strongly depend on the trapping geometry. 

Conclusion — We have established the boundary con- 
ditions at the interface between the SF core and the sur- 
rounding N shell for a trapped imbalanced fermion gas in 
an isotropic trap at ultralow temperatures. Using these 
conditions we have obtained the frequencies of collective 
modes, as a function of the polarization, when the normal 



gas behaves hydrodynamically. We have also calculated 
the frequencies of the monopole mode for a collisionless 
normal gas. 

For the hydrodynamic case, we find collective modes 
analogous to those in a single-component system, but 
with shifted frequencies. In addition, there exists a new 
class of collective modes; some of these have energies 
lower than the trap frequency. These modes correspond 
to out-of-phase motion of the SF core and the surround- 
ing N shell. 

For the case of a collisionless normal gas, we find that 
the lowest monopole mode remains at u> = 2wo, in agree- 
ment with an exact result due to Castin. We have also 
calculated the next lowest mode, finding that it reduces 
to the appropriate limits for vanishing and complete po- 
larization of the trap and, interestingly, that it is damped 
for some polarizations. This damping is a geometric ef- 
fect. It would be interesting to study crossover behavior 
of the OOP modes from the hydrodynamic to the colli- 
sionless normal gas. 
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